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Hydrogen bonding in infinite (HF)oo and (HCl)oo bent (zigzag) chains is studied using the ab 
initio coupled-cluster singles and doubles (CCSD) correlation method. The correlation contribution 
to the binding energy is decomposed in terms of nonadditive many-body interactions between the 
monomers in the chains, the so-called energy increments. Van der Waals constants for the two- 
body dispersion interaction between distant monomers in the infinite chains are extracted from this 
decomposition. They allow a partitioning of the correlation contribution to the binding energy into 
short- and long-range terms. This finding affords a significant reduction in the computational effort 
of ab initio calculations for solids as only the short-range part requires a sophisticated treatment 
whereas the long-range part can be summed immediately to infinite distances. 
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I. INTRODUCTION 

Hydrogen fluoride 12 ' 3 ' 4 ' 5 ' 6 and hydrogen chlorid o 7 ' 8 ' 9 
are representatives of molecular crystals; the electronic 
structure of the constituent HF or HC1 monomers is es- 
sentially preserved upon crystallization. The monomers 
in both crystals are hydrogen bondedfi°*ii*i2iiii4 it is a 
directional and anisotropic bonding of the hydrogen in a 
HF and HC1 monomer to the fluorine or chlorine atom, 
respectively, of a neighboring monomer. The bonding is 
caused by a partial withdrawal of charge from the hy- 
drogen atom due to the high electronegativity of the flu- 
orine or chlorine atoms liSjiLi^iiiiii Hydrogen bonds are 
intermediates between ionic bonding and van der Waals 
bonding and are of great importance for the physical and 
chemical properties of many organic and inorganic crys- 
tals. Moreover, they turn out to be crucial for the struc- 
ture of many biopolymers such as proteins and nucleic 
acids, 10 . 11 ! 12 . 13 ! 14 

At low temperature, HF and HC1 crystals are struc- 
turally very similar. While HF forms strong hydro- 
gen bonds, HC1 forms weak hydrogen bonds. There- 
fore, HF and HC1 represent good candidates for a thor- 
ough analysis of this special type of bonding in crys- 
tals. The monomers in both compounds are found to 
be arranged in terms of parallel zigzag chains with a 
large interchain distance and, hence, a weak interchain 
interaction [see Sec. 11111 for details]. Frequently, a sin- 
gle infinite chain is considered as a simple but realis- 
tic model of the crystals. The isolated (HF)^ chain 
has fascinated many theoreticians and much work has 
been carried out dealing with it; namely, model stud- 
iesfl£il& semiempirical (intermediate neglect of differen- 
tial overlap) examinations of (HF)^ by Zungerji 7 . den- 
sity functional theory calculations [local density ap- 
proximation] of Sprmgborg^* 1 ^ and ab initio investiga- 

tiOnSl 18 ! 19 ! 20 ! 21 ! 22 ! 23 ! 24 ! 25 ! 26 ! 27 ! 28 ! 29 ! 30 ! 31 ! 32 ! 33 In GOUteSt, 

the (HCl)oo chain has not been investigated very ex- 
tensively; there are only a few ab initio examina- 
tions i^liSii 3 ^ 5 . p or poth isolated chains the zigzag ge- 
ometry of Fig. ^ is the energetically favored arrange- 




FIG. 1: (Color online) Structure of infinite (HF) 

oo zigzag 

chains. The structure of (HCl)oo chains is analogous. The 
grey (red) numbers below the monomers indicate the rela- 
tive position of a monomer in the chain with respect to the 
monomer "0" at the origin. 



ment in comparison with the linear geometry. For (HF)oo 
chains this has been shown in early studies of Karpfen et 
qZ i 19 i 23 and Beyer and KarpfeniS 5 . The structure of iso- 
lated (HCl)oo chains is investigated in Ref. I3H 

The large number of ab initio studies that have been 
carried out for the chains also exhibits how challenging 
hydrogen-bonded systems are. An accurate treatment 
of hydrogen bonds requires both a good electron corre- 
lation method and a large one-particle basis set^S^ 
Hence, in Refs. I32l33l we applied basis set extrapolation 
schemes — which are powerful methods to improve the 
accuracy of both Hartree-Fock and correlation energies 
for small molecules — to extended systems, showing their 
validity and usefulness also in this situation. Thereby, 
a periodic Hartree-Fock treatment was accompanied by 
correlation calculations. 

We use the coupled-cluster singles and doubles (CCSD) 
correlation method^ 7 , in conjunction with the incremen- 
tal scheme^SiSMSiii^ which has been used in many 
applications and yields a physically meaningful many- 
body decomposition of the correlation energy. It has 
been applied to a few one-dimensional systems before: 
the ircms-polyacetylene polymer^ 3 " the infinite lithium 
hydride chain^ the beryllium hydride polymer^ and 
poly(para-phenylene)4 5 » In the context of the present pa- 
per, the contributions of Doll et a/.^^ 7 "* 4 ^ 9, who stud- 
ied ionic solids, as well as the investigations of rare-gas 
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crystals by Rosciszewski et alJ2l£i. are worth mentioning 
because hydrogen-bonded crystals fall in between these 
two types of bonding. 

This paper is structured as follows. The CCSD elec- 
tron correlation method and the incremental scheme are 
introduced in Sec. [H] whereas Sec. II I II describes geome- 
tries, basis sets, and the computer programs employed. 
In Sec. IIVI we discuss electron correlation effects and 
their implications. The contributions to the binding en- 
ergy of the (HFjoo and (HCl)oo chains are evaluated in 
Sec.lVland conclusions are drawn in Sec. I VII 



II. THEORY 

We choose the Hartree-Fock approximation as a start- 
ing point to study electron correlations, employing the 
full non-relativistic Hamiltonian of the chains within 
the fixed-nuclei approximationiSi 4 ^^ 2 **! 2 *^ 3 " The resulting 
Hartree-Fock Bloch orbitals are transformed to Wannier 
orbitals because the latter provide a more appropriate, 
local representation of the Hamiltonian for a subsequent 
treatment of electron correlations The ground- 

state wave function is modeled by the coupled-cluster 
ansataS 
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The coefficients i^ 1 " 1 and ^i ai >^2 a 2 arg re f erre( j ^ 
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as excitation amplitudes. Let r and s denote spatial 
and spin coordinates; then c 1 - creates electrons in 

virtual spin Wannier orbitals ws, (fs) whereas cs 

1 H x Ql ^ ' Hi Kl 

annihilates electrons from occupied spin Wannier or- 
bitals u>s (fs). Here R\ and R\ denote the unit cell 
in which the Wannier orbital is located and ot\ and ki 
refer to Wannier orbital indices. Assuming Born-von 
Karman boundary conditions, the correlation energy per 
unit cell £ corr is given by 



N S cmr = (*S\H\H!g)-($g\H\$$ 



(3) 



where Nq is the number of unit cells in the chains. 

Occupied Wannier orbitals are grouped in terms of n one 
pairwise disjunct one-body orbital sets which are defined 
by 



(1) 



which relates the iV-electron Hartree-Fock ground-state 
wave function $>q to the correlated ground-state wave 
function employing the cluster operator T. For the 
calculations in this study, it is sufficient to restrict T to 
single and double excitations, resulting in the coupled- 
cluster singles and doubles (CCSD) schem o 37 ! 55 



T 
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(2) 
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(4) 



for I = 1, 



The expression for the correlation energy of a crys- 
tal J2J is rearranged and decomposed in terms of cor- 
relation energies eg , s r of the electrons from the 

° Hi 1± - Rk Ik 

one-body orbital sets R\ 1% ■ ■ ■ Rk Ik ■ The resulting for- 
mula for the correlation energy of the chains per unit 
cell £ co „ reada 33 . 38 . 39 i 40 . 41 i 42 . 43 i 44 i 45 
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pairwise disjunct 

The one-body Ae^ t , the two-body Ae^ i r i ■> ■ ■ ■ > up to the i^-body Ae^ i 1 ...r k ■ ■ ■ energy increments are 
defined recursively fyr 33 , 38 ,^", 41 ! 42 , 43 ! 44 , 4 ^ 1 

Ae i?i h = £ Ri h ' ( 6a ) 
Ae Ri h 8.2 h = £ Ri h R2 I2 ~ Ae i?.i h " Ae R2 h ' ( 6b ) 
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The factors , A , . . . , , . . . in front of the sums on 
the right-hand side of Eq. JSJ account for permutations 
among the one-body orbital sets of a certain -K'-body en- 
ergy increment (K > 2). One can eliminate all permuta- 
tions that lead to the same energy increment by letting 
the sums in Eq. (J5J run only over distinct sets of one-body 
orbital sets. 

The translational relation of the Wannier orbitals leads 
to the translational symmetry of the energy increments, 
i.e., Aes r g s r ~ Aeg r g T . This can be ex- 

U 1\--Kk — H\ Ik Hi l\~-tlK Ik 

ploitcd in Eq. (JSJ to make the right-hand side of Eq. (J5J 
independent of the first lattice sum which, for this rea- 

son, is -/Vo times the sum of the translational symmetry- 
adapted energy increments. It allows us to eliminate the 
factor No in front of the left-hand side of Eq. JSJ) and enor- 
mously reduces the number of energy increments that 
need to be calculated to describe the chains with a given 
accuracy. 

The localized occupied molecular orbitals of oligomers, 
comprising a few monomers arranged in the geometry of 
the infinite chains, are found to be a good approxima- 
tion to the Wannier orbitals of infinite chains and crys- 
tals. 38 i 39 i 40 i 41 i 42 i 43 i 44 i 45 This facilitates the determination 
of the correlation contribution to the binding energy of 
infinite chains from the energy increments in oligomers 
using Eq. |JSJ. Yet the virtual Wannier orbitals involved 
in the determination of the energy increments, and im- 
plicitly contained in Eq. 10, are replaced by all the vir- 
tual canonical molecular orbitals of the oligomers. The 
procedure outlined in this paragraph is termed the incre- 
mental scheme, 38 i 39 i 40 i 41 i 42 



III. COMPUTATIONAL DETAILS 

Hydrogen fluoride and hydrogen chloride crystallize in 
an orthorhombic low-temperature phase described by the 
space groups Bm2 1 b for HF^L 2 - 3 . and Bb2 1 m for HClA 8 - 9 . 
The unit cells of both crystals contain four monomers 
which are arranged in terms of weakly interacting parallel 
zigzag chains [Fig. ^ . The chains are described by a unit 
cell which comprises two monomers and are considered as 
an excellent one-dimensional model for HF and HC1 crys- 

The structure of a single chain is determined by three pa- 
rameters, the H — X distance r, the X ■ ■ ■ X distance R, 
and the angle a = Z(HXH), X — F,C1. Experimental 
values for the parameters are r — 0.92 A, R = 2.50 A, and 
a = 120° for (HF)^^ and r = 1.25 A, R = 3.688 A, 
and a = 93.3° for (DCl)^. 7 HC1 and deuterated DC1 
crystals have very similar lattice constants and are 
considered to be isomorphousi Unfortunately, further 
structural information for HC1 crystals is unavailable. 



To calculate the energy increments in Eq. (J5J, we 
perform molecular calculations on oligomers (HF) n and 
(HC1)„, respectively, i.e., short fragments of the chains. 
We utilize the program package molpro, 5 - employing 
Foster-Boys localization^ 8 ^ of the molecular orbitals 
and the coupled-cluster singles and doubles (CCSD) cor- 
rclationmethodi2L££i§& Thereby, we are consistent with 
Refs. I32l33l where also the CCSD method was employed. 
We use the aug-cc-pVDZ basis se t 61 i 62 i 63 i 64 which yields 
somewhat less acc urate t otal binding energies than those 
obtained in Refs. I32I33L Yet we are mainly interested 
in the long-range behavior of electron correlations in the 
hydrogen-bonded chains; it is well described by the aug- 
cc-pVDZ basis set due to vanishingly small geometri- 
cal orbital overlaps among the one-body orbital sets in 
the -ftT-body energy increments considered. In fact, the 
modulus of the two-body energy increment involving the 
third- nearest-neighbor monomer |eo3|, as obtained with 
the aug-cc-pVDZ basis set, 61 i 62 i 63 i 64 is smaller by 2% 
for (HF)^ and by 4% for (HCl)oo than |e 3| obtained 
with the larger aug-cc-pVTZ basis set 61 i 62 i 63 i 64 [Sec. El- 
Similarly, if-body energy increments for K > 3 agree 
well as soon as the geometrical orbital overlaps between 
the one-body orbital sets become negligible. 



In Tab. IIIII below, reference is made to periodic 
Hartree-Fock calcu lations of (HF)oo and (HCl)co chains, 
described in Refs. 32 33. which were carried out with 
the crystal program. 65,66 '— For these calculations, the 
f,g,h,i functions were removed from the basis sets. To 
estimate the influence of the neglected basis functions 
on the total binding energy of the chains, we use MOL- 
PRO)££ to calculate the Hartree-Fock binding energy per 
monomer for (HF)g and (HCl)g both with and without 
/, <?, h functions using the cc-pV5Z basis setj^LSiSi at the 
Hartree-Fock level it is close to completeness with respect 
to basis functions with the angular momenta s,p, d. The 
binding energies of both chains obtained with f,g,h func- 
tions are smaller by 13 /iEh than those determined 
without them. This deviation turns out to be sufficiently 
small to be neglected in the subsequent calculation of to- 
tal binding energies in Sec. However, note that for 
the correlation calculations the /, g, h functions are fully 
considered. 



The calculated binding energies reported here account 
for the basis set superposition error— (BSSE) by the 
counterpoise correction (CP); it is determined by sur- 
rounding an isolated HF or HC1 monomer by sufficiently 
many ghost atoms positioned in the geometry of the in- 
finite chains&g 9 . 
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FIG. 2: (Color online) One-body energy increments Aeo 
of (HF)oo and (HCl)oo chains as determined in oligomers of 
varying length reduced by the CP-corrected correlation en- 
ergy £ mon °f the corresponding isolated monomer. The mean 
value of the two innermost energy increments of the oligomers 
is taken. Circles stand for (HF)oo and squares for (HCl)oc 
data. Open symbols refer to oligomers surrounded by point 
charges whereas closed symbols denote isolated oligomers. 



IV. ELECTRON CORRELATIONS 
A. Transferability 

We have to ensure that the Foster-Boys-localized 
molecular orbital o 58 i 59 which are extracted from 
oligomers are a good approximation to the Wannier or- 
bitals of the infinite (HF)^ and (HCl)oo chains. They 
have to be approximately translationally related within a 
certain oligomer and must not differ substantially among 
oligomers of varying length. These two properties of the 
localized occupied molecular orbitals are termed trans- 
ferability. In order to achieve transferability for moder- 
ately sized molecular clusters, one is frequently obliged to 
account for the omitted monomers in terms of an appro- 
priate substitute. In the case of (HF)^ and (HCl)oo, two 
options offer themselves. On the one hand, both chains 
are cut out from a molecular crystal, implying no sub- 
stitutes as in rare-gas crystals on the other hand, 
particularly (HF)^ is rather ionic and a surrounding by 
point charges can be envisaged to create the embedding 
which adequately models the infinite chains4S^£^2^2iBi 
I explore both possibilities to elucidate what kind of pro- 
cedure is adequate for hydrogen-bonded solids. The vir- 
tual orbitals of the oligomers are not altered and remain 
canonical molecular orbitals in the calculation of the indi- 
vidual energy increments. The oligomer approximation 
thus corresponds to a sort of domain decomposition of 
the virtual space which also is applied in other local cor- 
relation methods^ 

In the case of (HF)^ and (HCl)^, the outer and in- 
ner valence orbitals, i.e., eight electrons, are assigned to 



FIG. 3: (Color online) Two-body energy increments Aeoi 
of (HF)oo and (HCl)oo chains as determined in oligomers of 
varying length. The two innermost one-body orbital sets of 
the respective oligomer are taken to determine Aeoi • Symbols 
are chosen as in Fig. |5] 



the one-body orbital set representing a particular HF or 
HC1 monomer (0J). To test the transferability, we calcu- 
late all one-body energy increments and all two-body en- 
ergy increments between two adjacent monomers, the so- 
called connected two-body energy increments, in isolated 
oligomers (HF)„ and (HC1)„ for n = 2, 4, 6, 10. They are 
compared with the energy increments from the oligomers 
which were surrounded by point charges up to twentieth- 
nearest neighbors. The hydrogen atoms are represented 
by q = +1; fluorine and chlorine atoms are described 
by q — — 1 with the exception of the terminal charges 
which are set to q = - 1,46,47,48,49,70 

Comparing the one-body and the connected two- 
body energy increments obtained in oligomers of varying 
length in Figs. and we observe a rapid convergence 
of their values toward the limit of the infinite chains 
where the values determined for the isolated oligomers 
and the values resulting from the point-charge-embedded 
oligomers approach each other quickly. The energy in- 
crements taken from the isolated oligomers increase (de- 
crease) in Fig. |21 (Fig. |2J monotonically with the oligomer 
length whereas the embedding with point charges causes 
an unsystematic but quicker convergence. The conver- 
gence behavior of the individual curves in Figs. [5] and 
is governed by a small variation in the localized orbitals 
and, predominantly, by the improvement of the electronic 
structure which approaches rapidly the electronic struc- 
ture of the infinite chains with increasing length of the 
oligomers. Obviously, the decision whether to employ 
point charges or not has only a small impact on the val- 
ues of the energy increments. The one-body (two-body) 
energy increments differ by 67 fiEh (14/zE/j) in (HF)^ 
and by 22 fiE^ (31 fiEh) m (HCl)oo. In what follows only 
energy increments taken f rom isolated oligomers are con- 
sidered, as in Refs. I32I33L 
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Increment (HF)oo (HCl)oo 

As~ -223.3632 -175.0985 

Ae i -4.0879 -2.8213 

Ae 02 -0.0861 -0.1347 

Ae i2 -0.0163 0.0014 

Ae 03 -0.0063 -0.0092 

Ae i3 -0.0015 -0.0016 

Ae i23 0.0003 0.0006 



TABLE I: Exemplary energy increments of (HF)oo and 
(HCl)oo chains taken from (HF)io and (HCl)io oligomers. The 
energy increments are computed using the one-body orbital 
sets of the innermost monomers in the oligomers. All data 
are given in millihartrees. 



The error on the binding energy due to the approxi- 
mate translational relation of the Wannier orbitals taken 
from the finite-length oligomers can be estimated by in- 
vestigating the numerical differences between energy in- 
crements whose values should be identical due to (trans- 
lational) symmetry. The connected two-body incre- 
ments, for example, vary at most by 40 /iEh in (HF)io 
and at most by 11 /iEh in (HCl)io^ We conclude that 
the transferability prerequisite of the energy increments 
is satisfied within good accuracy. The agreement be- 
tween the results obtained in the two different ways de- 
scribed before corroborates the applicability of oligomers 
to represent the electronic structure of infinite (HF)oo 
and (HCl)oo chains. 



B. Short-range correlations 

Let us discuss the short-range correlation contributions 
first. The most relevant energy increments of (HF)^ and 
(HCl)oo are summarized in Tab. [I] The CP-corrected cor- 
relation energy of the HF monomer, —227.3406 mE^, and 
of the HC1 monomer, —176.2151 mE^, are only larger in 
magnitude by 1.8% and 0.6%, respectively, than the cor- 
responding one-body energy increments Asq. They are 
a bit larger because the electrons of a monomer in the 
chains experience Pauli repulsion exerted by the neigh- 
boring monomers. However, these effects are small in 
our case and the electronic structure of the HF and 
HC1 monomers is essentially preserved in the (HF)oo and 
(HCl)oo chains. 

The -ftT-body energy increments, K > 2, describe the 
mutual correlation of the valence electrons of several 
monomers. They give rise to a pronounced nonlinear 
increase of the binding energy of small clusters and short 
oligomers which is termed bond cooperativity (Ref. ^| 
and references therein). The modulus of the connected 
two-body energy increment of (HCl)oo, |Ae^° l |, is 31% 
smaller than the modulus of the corresponding energy 
increment of (HF)oo, |A£qj |. The connected three-body 
energy increment Ae^j 1 is repulsive, i.e., greater than 
zero, but Ae^ is again attractive, i.e., smaller than 



Compound Cg dd Cg vcn Experiment" 



(HF)oc 


23.4 


20.2 


19.00 


(HCl)oo 


121 


120 


130.4 



"The experimental data are taken from Ref. [7S it is a rotational 
average of the van der Waals constants for the dispersion interac- 
tion between two isolated monomers. 



TABLE II: Van der Waals constants and C£ ven of the 

two-body dispersion interaction for the (HF)oo and (HCl)oo 
chains. They are obtained by fitting the two-body energy in- 
crements of Fig.0]with either odd or even cardinal numbers n, 
respectively, to Eq. 0. All data are given in hartree bohr B . 



zero. The reverse trend is observed for the remaining 
energy increments. The modulus of the energy incre- 
ment Ae[/ 2 C1 is 57% larger than |Ae^ 2 F |. The connected 
four-body increment Ae^g is 100% larger than Ae^ 3 
and both are repulsive. These two trends of the en- 
ergy increments can be explained by two effects. First, 
short-range correlation is effective for nearest neighbors 
[see also the ensuing Sec. II V Cj : it is apparently stronger 
in (HF)oo than in (HCl)oo due to the tighter packing of 
the monomers in (HF)oo and the higher compactness of 
the HF monomer itself. Second, chlorine atoms have a 
higher polarizability than fluorine atoms because their 
valence electrons are more diffuse than those of the lat- 
ter atom. This causes the van der Waals interaction to 
be stronger in (HCl)oo than in (HF)oojS leading to the 
more distant energy increments being larger in (HCl),^ 
when compared with (HF)^; all this despite the larger 
intcrmonomer distances in (HCl)oo than in (HF)^. 

However, note that one should not associate too much 
physics with the particular values of energy increments; 
namely, they depend on the specific unitary transforma- 
tion used to localize the occupied orbitals (here the one of 
Foster and Boys^S) . Only the binding energy, which in- 
volves the sum of all the energy increments [Eqs. © and 
(|10f) ]. is a physical observable and thus invariant under 
orbital transformations. When an analysis of the indi- 
vidual energy increments is to be meaningful, their value 
should be fairly independent of the localization procedure 
which is well satisfied in our case. 



C. Long-range correlations 

At a separation of two isolated monomers where the ge- 
ometrical overlap between the orbitals from distinct one- 
body orbital sets of the monomers becomes negligible, 
only van der Waals dispersion interactions remain. This 
is also the case for the interaction between two monomers 
in the infinite chains which is correspondingly described 
by the two-body energy increments. Their absolute val- 
ues in (HF)oo and (HCl)oo for two monomers at a dis- 
tance up to seventh nearest neighbors are displayed in 
Fig. the curves drop off rapidly with distance. 

The monomers in Fig. ^ that are labeled by odd cardi- 
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nal numbers are tilted with respect to monomer 0. Like- 
wise, monomers with even cardinal numbers are arranged 
parallel to the monomer 0. We will refer to the two types 
of two-body energy increments that result from a parallel 
or tilted setting of monomers as even or odd energy incre- 
ments, respectively. The long-range interaction between 
two monomers is approximated by the leading term of 
the two-body van der Waals dispersion interaction 72,74 



1.0x10 



_vdW 
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■ifw ' n even ' 



(7) 



where individual van der Waals constants Cg and 
C| ven , respectively, are affixed for odd and even energy 
increments. Here a represents the lattice constant of the 
chains. The van der Waals constants are obtained by a 
weighted nonlinear curve frfeS of the data in Fig. 0] for ei- 
ther the odd or the even energy increments. The weights 
are chosen such that the fit reproduces the energy in- 
crements with large n best, as for them the geometrical 
orbital overlaps between the orbitals from two different 
one-body orbital sets are negligible. The van der Waals 
constants for the fits of the two sets of translationally 
equivalent monomers are given in Tab. [H] Our theo- 
retical data compare satisfactorily with the experimental 
data for the dispersion interaction between two isolated 
monomers also given in the table. 

Having identified the van der Waals contribution to the 
two-body energy increments in the chains, we can sub- 
tract it from the two-body energy increments to obtain 
van der Waals-reduced energy increments 



n.e . 
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(8) 



Their absolute value | Aeq„ | is shown next to the abso- 
lute value of the two-body energy increments |Aeon| in 
Fig. m 

Two regions can be identified in Fig. ^ First, from 
the nearest to the third-nearest neighbors, there is the 
local correlation zone, where, due to geometrical orbital 
overlaps, short-range electron correlations are effective. 
There, the decay of two-body energy increments with the 
distance between the two monomers is faster than the one 
that would result from a pure van der Waals interaction, 
following Eq. 10. We observe that the van der Waals 
contribution to A£p dw | of (HCl)oo is appreciably larger 
than in (HF)^ as the curves for |Ae i| and |AeQ dw | are 
much closer for (HCl)oo than for (HF),^. Second, there 
is the van der Waals zone reaching from the third- up to 
the seventh-nearest neighbor. Here a typical r decay is 
observed, leading to van der Waals-reduced energy incre- 
ments that are essentially zero. A slight deviation of the 
two-body energy increments from an r~ 6 behavior is per- 
ceived beyond fifth-nearest neighbors which likely can be 
attributed to inaccuracies caused by the oligomer approx- 
imation. The translational relation of the Wannier or- 
bitals is better satisfied in (HCl)oo compared with (HF)oo 
[Sec. IIV A| which leads to a lower absolute value of the 




FIG. 4: (Color online) Absolute values of two-body energy in- 
crements |Aeon| (closed symbols) and van der Waals-reduced 
two-body energy increments |A£gn W | (open symbols) using 
Eq. and the constants in Tab.lTTlfor both the (HF)oo chain 
(circles) and the (HCl)oo chain (squares). 



fourth to seventh van der Waals-reduced energy incre- 
ments. The low absolute value of AeQ dw for both chains 
is an artifact of the fitting process. Figure 0] also reveals 
that the two-body energy increments are of satisfactory 
accuracy even beyond the estimates given in Sec. IIV Al 
of 40 ^E h for (HF)io and 11 ^E h for (HCl)i due to a 
considerable error cancellation. 

Given the van der Waals constants in Tab. |H] we 
can use Eq. to sum up the van der Waals contri- 
bution of the two-body energy increments to the bind- 

oo 

ing energy to infinite distance ^ ^on W which yields 

n— 3 

-8/xEh for (HF)^ and -12 ^E h for (HCl)^. The bigger 
polarizability of chlorine atoms compared with fluorine 
atoms leads to the larger contribution in (HCl)oo when 
compared with (HF)^ as can be expected from the anal- 
ysis of the energy increments in Sec. IIV 51 The contribu- 
tion of two-body energy increments beyond third-nearest 
neighbors is in both cases ~ 30% of Aeo3- 

Generally, once a good estimate of the van der Waals 
constants of a crystal is available, the decomposition of 
the energy increments into a short- and a long-range van 
der Waals contribution can be made. This allows one 
to focus on the short-range part. Thereby, the required 
number of energy increments to be calculated to reach 
a certain accuracy of the correlation energy is reduced 
considerably. Therefore, the computational demand of 
its determination is also significantly decreased. 



V. BINDING ENERGY 

Having understood the impact of electron correlations 
in (HF)oo and (HCl)oo chains, we are now in a position 
to examine the total binding energy per monomer AE of 
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the chains; it is given by 
AE = AE SCF + AE COII 

_ 1 / cSCF i ccorr \ _ cSCF _ ccorr 

2 v chain * chain/ ^monomer ^monomer ' 

(9) 

It consists of the Hartree-Fock and the correlation con- 
tributions ASgcF and AE COII , respectively, where ffhafn 
and £chain denote the Hartree-Fock and correlation ener- 
gies per unit cell of the chains. Furthermore, Monomer 
and ^monomer are the corresponding energies for the 
monomers. The Hartree-Fock energies are readily avail- 
able from molecular and periodic calculationsSiSi*§! but 
electron correlations are more involved and were dis- 
cussed in the previous Sec. IIVI 

In Tab. IIIII we communicate basis-set-extrapolated 
Hartree-Fock binding energies A_Escf(oo) for the chains 
taken from Refs. At the Hartree-Fock equi- 

librium geometry of (HF)oo, Bayer and Karpfen 
give — 10.360mE/i for basis set 3 in Ref. I25I while Hirata 
and Iwata find -10.855 mE h for the 6-311++G(d,p) ba- 
sis set^£ Both numbers are in good agreement with our 
result which is — 10.202mE/,,. Yet Berski and Latajka re- 
port —9.696 mE/j for (HF)^, when using the 6-311G(d,p) 
basis set^ and — 2.073 mE^ for (HCl)oo harnessing the 
DZ(d,p) basis setm^ The latter number in particular de- 
viates considerably from our result, —0.586 mE/,, which 
is most likely due to the fact that they did not remove 
the BSSE,^ 

To determine the correlation energy per monomer of 
the infinite chains £ C orr, we inspect the energy increments 
in Tab.[I] They converge rapidly both with increasing dis- 
tance among the monomers involved and with increasing 
order of the many-body expansion, i.e., with the max- 
imum number of monomers correlated in a specific en- 
ergy increment. In order to obtain an accuracy in £ corr 
of fa 50/zE h for (HF)^ and of « 11 /xE h for (HCl)^, we 
see from Tab.Uthat it is enough to include the following 
energy increments: 

£ corr ra Ae + Ae i + Ae 02 . (10) 

Hence, it is sufficient to correlate the electrons in 
the Wannier orbitals of only two monomers at a 
time for (HF)^ and (HCl)oo chains in which more 
distant two-body energy increments contribute only 
— 8 fiE h for (HF)^ and — 12^iE h for (HCl)^ as shown 
in Sec. IIV CI In three-dimensional crystals, three-body 
terms become more important according to the analyses 
111 Refs. 151171771 and the much larger number of three- 
body energy increments compared with the number of 
two-body energy increments. 

We report AE based on the approximation of Eq. i|10|) 
as obtained by basis-set extrapolation^^ in Tab. IIIII 
Using density functional theory, Hirata and Iwata™ 
find at optimized geometries for the total binding en- 
ergy of (HF) M -13.521 mE h (BLYP) and - 13.864 mE^ 
(B3LYP) [utilizing the 6-311++G(d,p) basis set] which 
are larger in magnitude by 14% and 17%, respectively, 
than our result for the experimental geometry. The ob- 





(HF)oc 


(HCl)oo 


A_E S cf(oo) 


-10.202 


-0.586 


A-Ecorr(oo) 


-1.623 


-2.764 


AE(oo) 


-11.826 


-3.350 



TABLE III: Basis-set-extrapolated total binding energies 
per monomer AE(oo) of (HF)oo and (HCl)oo chains and 
their decomposition into basis-set-extrapolated Hartree- 
Fock A_Bscf(oo) and electron correlat ion A _E corr (oo) contri- 
butions. All data are taken from Refs. I32i33l and are given in 
millihartrees. 



served differences can be partially ascribed to the differ- 
ent geometries used; our experimental geometry leads to 
a binding energy that is somewhat smaller in magnitude 
than the binding energy in the optimized geometry. The 
6-311+- |-G(d,p) basis set of Iwata and Hir ata^ a lso dif- 
fers substantially from our choices in Refs. l32l33l Above 
all, we do not consider the impact of triple excitations in 
our study; they increase the absolute value of the short- 
range correlation contributions and, consequently, also 
the absolute value of the binding energies. 



VI. CONCLUSION 

We study hydrogen bonding in infinite (HF)oo and 
(HCl)oo chains. Whereas in (HF)oo the Hartree-Fock 
contribution dominates the total binding energy by 86% 
due to the electrostatic contribution of the rather ionic 
HF monomers, in the much more weakly bound (HCl)oo, 
the Hartree-Fock calculation yields only 18% of the bind- 
ing energy [Tab. IHIj . The transitional character of 
hydrogen-bonded crystals between the ionic and van der 
Waals regimes is reflected in these numbers. The very 
weak hydrogen bonds in the (HCl)oo chains bear a close 
resemblance to purely van der Waals-bonded systems, 
like rare-gas solids, where bonding is entirely caused 
by electron correlations In fact only the inclu- 
sion of electron correlations puts the binding energy per 
monomer of the (HCl)oo chain into the range convention- 
ally ascribed to hydrogen bonding^! namely 3-16 mE^. 

The incremental scheme, which provides a decom- 
position of the contribution of electron correlations to 
the binding energy of the chains in terms of nonad- 
ditive many-body energy increments, is shown to con- 
verge rapidly with respect to the number of monomers, 
i.e., bodies, involved in the energy increments and the 
distance among them, thus providing a good tool to 
study hydrogen bonding in crystals. In contrast to 
three-dimensional crystals, the dominant contribution 
of electron correlations to the binding energies of the 
infinite chains is, in the present investigation, already 
given by pair interactions between the valence electrons 
of a monomer and its nearest and next-nearest neigh- 
bors. The if -body energy increments represent a gen- 
eral framework to describe van der Waals dispersion in- 
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teraction among the electrons from K one-body orbital 
sets. The decay of the two-body energy increments, i.e., 
the electron correlations between two monomers in the 
chains, with the intermonomer distance is investigated 
and van der Waals constants are determined, which af- 
fords a partition of the energy increments into a short- 
range, van der Waals-reduced part and a long-range part 
that can be summed analytically to infinite distances. 

This partitioning affords a methodological advance- 
ment; if the van der Waals constants of a crystal are 
available from other sources, then the number of energy 
increments in the incremental expansion can be reduced 
by considering only van der Waals-reduced energy incre- 
ments; they typically need to comprise only the nearest- 
and next-nearest-neighbor energy increments. The small 
number of van der Waals-reduced energy increments per- 
mits one to significantly decrease the number of energy 



increments which have actually to be calculated to reach 
a certain accuracy of the correlation energy compared 
with the conventional incremental series. Furthermore, 
the long-range part turns out to be well described us- 
ing small (double-^) basis sets — in contrast to the short- 
range part where basis-set extrapolation is required for an 
accurate descriptions^^ — thus facilitating further com- 
putational savings. 
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